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Abstract. Recently Dantas, Oliveira and Stilck [J. Stat. Mech. (2007) P08009] 
studied how the one-dimensional diffusive contact process crosses over from the critical 
behavior of directed percolation to an effective mean field behaviour when the diffusion 
rate is sent to infinity. They showed that this crossover can be described in terms of 
a crossover exponent finding the boundaries 3 < (/) < 4 in one spatial dimension. In 
the present work we refine and extend this result up to four spatial dimensions by a 
field-theoretic calculation and extensive numerical simulations. 



PACS numbers: 05.50.+q, 05.70.Ln, 64.60.Ht 
1. Introduction 

In non-equilibrium statistical mechanics, the study of phase transition from fluctuating 
phases into absorbing states continues to be a very active field [1]. One of the most 
important universality classes of such transitions is directed percolation (DP) [2-5], 
which is well understood and can be described in terms of a renormalizable field 
theory [6,7]. Although this universality class plays a paradigmatic role as the Ising 
class in equilibrium statistical mechanics, it was primarily of theoretical interest since 
experimental realizations seemed to be very difficult [8]. However, very recently Takeuchi 
et al. [9] performed an experiment on the basis of turbulent liquid crystals, where the 
critical exponents of directed percolation could be measured. 

Directed percolation is most easily introduced as a spreading process of diffusing 
particles on a lattice which multiply {A 2 A) and self-annihilate (A — 0). Moreover, 
the density of particles is effectively limited by an exclusion principle. One of the most 
important models with random-sequential updates is the contact process [1] which is 
controlled by the rate for offspring production A with a critical point Ac. The DP 
universality class is characterized by three critical exponents z/y, i/^) which depend 
only on the dimensionality of the system. In = 1, 2, 3 spatial dimensions the exponents 
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take non-trivial, probably irrational values. In > 4 dimensions, however, the mean- 
field exponents /3 = z/y = 1 and v± = 1/2 become exact. Right at the upper critical 
dimension dc = 4: there are additional logarithmic corrections. 

The crossover to mean field behavior for d > 4 is related to the fact that the diffuse 
mixing becomes more efficient in high dimensions, suppressing fiuctuation effects. An 
alternative way to enhance diffusive mixing, on which we will focus in the present work, 
is to increase the diffusion constant D in a low-dimensional system. In most models 
such as the contact process, diffusion is effectively included by the circumstance that 
a particle creates offspring at an empty randomly chosen nearest-neighbour site. In 
such models the effective diffusion rate is fixed but it is straight forward to add explicit 
diffusion of solitary particles so that the diffusion rate can be controlled. 

By varying the diffusion constant in the contact process, one expects the following 
phenomenology. Clearly, in the limit D ^ oo, where all sites are mutually coupled, mean 
field theory becomes exact. For a large but finite diffusion constant, however, one expects 
a crossover at a typical length scale ^^f^ and an associated typical time scale ^^^^ which 
grow with D. Below these scales, the process exhibts an effective mean-field behavior 

(c) (c) 

while it crosses over to DP on scales larger than and ^ . In addition, the critical rate 
for offspring production, Ac(-D), decreases with D and reaches the analytically known 
mean field value Ac(cxo) in the limit D — > oo. 

The infiuence of the diffusion rate on Ac was first studied by Konno [10] from a 
mathematical point of view. He pointed out that the critical threshold Ac(-D) approaches 
the asymptotic value Xc{oo) algebraically as 

A,(D) - Ac(oo) ~ D-'/"^ , (1) 

where the crossover exponent </> = 3 in one, = 1+ log. corrections in two and = 1 
in more than two spatial dimensions. More recently, by series expansion and partial 
differential approximants, Dantas et al. [11] found the boundaries 3 < < 4 for this 
exponent. In the present work we describe a pure field theoretic approach which allows 
an easy calculation of the crossover exponents. Furthermore we present numerical results 
for one to four spatial dimensions. 



2. Field-theoretical approach 

The diffusive DP model is most easily described by a contact process defined by the 
three micro processes (i) creation of a particle on, (ii) hopping to one of the neighboring 
places, and (iii) death of a particle which take place with rates A, D and 1 respectively. 
These dynamic laws result in the coarse grained Langevin equation [12] 

|-p (f , t) = ^ Ap (x, t) + (A - 1) p (x, t) - Ap2 (f, t)+i (f , t) (2) 
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with spatial dimension d, density of particles p {x, t) and multiplicative noise ^ (x, t) 
with the correlations 

(e (x, t) i {y, t')) ~ p (f , t) 5" {x-y)5{t- t') . (3) 

Based on this equation one can set up a field theory by eliminating the noise and 
introducing a response field. After appropriate rescaling of fields and parameters, the 
field theoretic action reads (see [7] and references therein) 

S = j d'xdt^ (^0^0 - 0^0^ +^ (V0) (V0) -4>f^<P+g - 4>'<p) (4) 

with density (response) fields (0), control parameter k = A — 1 and coupling g. 
As seen from ([4]) the field theory of DP involves a directed propagator and two cubic 
interaction vertexes with the same coupling strength g. Using these Feynman rules, 
one can calculate the next order (one loop) corrections to propagator and vertexes. For 
instance, the inverse propagator to one loop order is given by 

T^'^'Hk,uj)= H'^.- +0{g'). (5) 

Naturally these loop corrections lead to divergent contributions, which receive a physical 
meaning upon applying a regularization and renormalization procedure. The divergences 
then appear as corrections to the meanfield scaling exponents 

^± = ^ + ^ + (^(e)% u\\ = l + ^ + 0{e)\ f]=l-'- + Oie)' (6) 

where e = 4 — d. As expected, the critical exponents do not depend on the value of the 
diffusion rate. 

To get an insight how the diffusion rate D changes the critical creation rate Ac, let 
us consider the field theoretic calculation at the beginning. Due to convergence problems 
upon applying dimensional regularization, one is forced to replace the parameter n by 
mass = k — (D) where 

f^c = g^ [ with ^j2dKjD <\q\<n (7) 

is chosen in a way that m becomes zero at criticallity [7]. The lower boundary is given 
by the fact, that for q smaller than ^J2d^^c/ D the propagator becomes zero. The upper 
boundary f2 is a cutoff scale in momentum space and should be sent to infinity under 
continued renormalization. Inserting the rf-dimensional surface element and expanding 
the right hand side of (JT]) in powers of q one obtains an integral over a geometric series 



Kr. 



DT {d/2) J^y^;^ V ^ 
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Figure 1. Log-Log plot of Xc{D) derived by Monte-Carlo simulation. The straight 
lines show power laws, fitted to the large D part of the numerical data. 



By dimensional analysis one can show that for d < 2 only the lower boundary and for 
d> 2 only the upper boundary contributes to the integral. As we are interested in the 
asymptotic behavior for large we may approximate the sum by the leading term 1. 
Integrating the remaining part then yields the crossover behavior for different spatial 
dimension 

{D^^d for d < 2 
M for = 2 (9) 
D-^ for d > 2 

These results compare very well to the predictions made in [10] but here they are 
obtained in a much simpler and direct way. From Q we read off the crossover exponents 

{3 for c/ = 1 

1 + log. corrections for = 2 (10) 
1 for c/ > 3 



3. Numerical results 



In order to substantiate these analytical results, we simulated the diffusive contact 
process on periodic rf-dimensional lattices using the following dynamic rules: Select a 
random lattice site and if there is a particle perform one of the following moves: 

(i) remove the selected particle with probability 1/ (1 + D + A), 



Crossover from directed percolation to mean field behavior 5 

10 

1 

0.1 
0.01 
0.001 



1e-04 

0.1 1 10 100 1000 

D 

Figure 2. The same simulation results as in figure [1] compared to the numerically 
integrated predictions of eq. ([7]) (solid lines). Since the iteration becomes unstable for 
small diffusion rate, this region is not accessible by this approach. 




(ii) create a particle on a randomly chosen neighbouring site with prob. A/ {1 + D + \) 
if this site is empty, or 

(iii) move the particle to randomly selected neighbouring site with prob. D/ {1 + D + \) 
if this site is empty. 

For various values of the diffusion rate D we measured the particle density in a system 
starting with a fully occupied lattice and averaged over at least 20 independent runs. 
By varying the creation rate A and searching for a power law decay we determined the 
critical threshold Ac as a function of D. These critical lines for one, two, three and four 
spatial dimensions are shown in figure [H The plots clearly indicate the expected power 
law behavior for large diffusion rates and the exponents are in good agreement with our 
analytical results. For small D, however, the observed deviations from a pure power law 
are not very surprising because the diffusion length becomes comparable to the lattice 
spacing. 

So far we verified the predicted asymptotic behavior ([9]). In order to investigate 
the next-leading corrections of the geometric series, we integrate ([HI) numerically and 
compare it with the results of the simulations (figure [2]). Although for small D this 
solution is plagued by numerical instabilities rendering the results unusable, the reader 
should notice that for spatial dimension d = 1 the numerical integration becomes almost 
exact over at least two decades and does not show the lift off behavior like in other 
dimensions. 
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4. Concluding remarks 

In this paper we have presented a simple field-theoretic calculation and extensive 
numerical simulations in order to investigate how the critical parameter depend on the 
diffusion rate of a DP process in various dimensions. Our field-theoretical result confirms 
previous results by Konno [10], derived here in a much simpler way. Moreover we 
have presented numerical results which are in agreement with these predictions, refining 
previous results by Dantas et al. [11] and extending them to higher space dimensions 
up to (i = 4. Especially in one spatial dimension, the field-theoretic prediction to one 
loop order, when evaluated numerically, coincides almost perfectly with the Monte Carlo 
estimates. This is surprising since one expects loop corrections to be more relevant in 
lower dimensions. 
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